Pion induced double charge exchange reactions in the A resonance region 
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We have applied the Giessen BUU (GiBUU) transport model to the description of the double 
charge exchange (DCX) reaction of pions with different nuclear targets at incident kinetic energies 
of 120 — 180 MeV. The DCX process is highly sensitive to details of the interactions of pions with 
the nuclear medium and, therefore, represents a major benchmark for any model of pion scattering 
off nuclei at low and intermediate energies. The impact of surface effects, such as the neutron skins 
of heavy nuclei, is investigated. The dependence of the total cross section on the nuclear mass 
number is also discussed. We achieve a good quantitative agreement with the extensive data set 
obtained at LAMPF. Furthermore, we compare the solutions of the transport equations obtained in 
the test-particle ansatz using two different schemes - the full and the parallel ensemble method. 
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I. INTRODUCTION AND MOTIVATION 

Over the last two decades, the Giessen Boltzmann-Uehling-Uhlcnbcck (GiBUU) transport model has been developed 
to describe heavy ion collisions, photon-, electron-, pion- and neutrino- induced reactions within one unified transport 



framework. A realistic treatment of pion interactions with the nuclear medium is crucial for the interpretation 
of experiments where pions are produced inside nuclei. Recently, it has been shown using GiBUU 0, S i that 
pion rescattering in the final state description of photon induced double-pion productionproduce a considerable 



pion rescattering m tne nnai state description oi pnoton induced double-pion production produce a considerable 
modification of the 7T7t invariant mass distributions observed by the TAPS collaboration Moreover, neutrino- 
induced pion production, a source of background for neutrino oscillation experiments is very sensitive to pion final 
\q \ state interactions |6J, |j] . The propagation of low-energy pions in nuclear matter in the GiBUU framework has already 
■ been extensively discussed in ref. [2| and compared to quantum mechanical calculations. 

In this context, pionic double charge exchange (DCX) is a very interesting reaction. The fact that DCX requires 
at least two nucleons to take place makes it a very sensitive benchmark for pion rescattering and absorption. This 
reaction received a considerable attention in the past (see for instance Ref. Q and references therein). The mechanism 
of two sequential single charge exchanges has traditionally been able to explain the main features of this reaction 
at low energies although the contribution of the A(ir, irn)X reaction becomes progressively important as the energy 
increases 

[ill O At hi g her L± 1 GcV ) 

energies, the sequential mechanism becomes insufficient to account for 
the reaction cross section 0, ^| . Extensive experimental studies performed at LAMPF obtained high precision 
data for doubly differential cross sections on 3 He [13 and heavier nuclei ( le O, 40 Ca, 208 Pb) 16] in the region of 
E kin = 120 - 270 MeV. 

Hiifner and Thies |l7j explored for the first time the applicability of the Boltzmann equation in ttN collisions and 
achieved qualitative agreement with data on single and double charge exchange. Their method to solve the Boltzmann 
equation was based upon an expansion of the pion one-body distribution function in the number of collisions. There, 
in contrast to our work, the Boltzmann equation is not solved with a test-particle ansatz but by reformulating it into 
a set of coupled differential equations which can then be solved in an iterative manner. However, this approach was 
based on simplifying assumptions of averaged cross sections and averaged potentials. The work by Vicente et al. [Tlj 
was based upon the cascade model described in [l^. There, a microscopic model for nN scattering was used as input 
for the pion reaction rates in the simulation. In that work pion DCX off 16 and 40 Ca has been explored and 
fair quantitative agreement with data has been achieved. 

In our work, we explore DCX on heavier nuclei, comparing with the data measured by Wood et al. ^(|. We also 
address the scaling of the total cross section discussed by Gram ct al. . To focus only on single-pion rescattering, 
we consider incoming pion energies below Skin = 180 MeV; above that energy 27r production becomes prominent and 
DCX does not happen necessarily in a two-step process anymore. Due to the small mean free path of the incoming 
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pions, the process is mostly sensitive to the surface of the nucleus. Therefore, we will discuss and compare two 
widely used numerical schemes for the solution of the Boltzmann equation: parallel ensemble method employed in the 
BUU models Etl l2ll I221 123^ and in the Vlasov-Uehlin g-U hlenbeck model [24|; and full ensemble method used in the 
Landau- Vlasov[25l| . Boltzmann- Nordheim-Vlasov [26l |27| and Relativistic BUU [23, |2^| models. Both schemes are 
based on the test-particle representation of the single-particle phase space density, but they differ in the locality of 
the scattering processes (c.f. discussion in |30l l3lj ). In the parallel ensemble method, all test-particles are subdivided 
into the groups, or parallel ensembles. The number of test-particles in each parallel ensemble is equal to the number 
of physical particles in the system. Collisions are allowed only between test-particles from the same parallel ensemble, 
while the mean field is averaged over all parallel ensembles. Without mean field, thus, the parallel ensemble method 
is equivalent to the intranuclear-cascade simulation |32l |. In the full ensemble method, collisions between all test- 
particles are allowed. The low-energy (E^ = 10—50 MeV/nucleon) heavy ion collisions are better described by the 
full ensemble method, since this method provides convergence to the exact solution of the original kinetic equation in 
the limit of large N = (number of test-particles per nucleon). However, in codes oriented to the particle production 
in high-energy heavy ion and hadron-nucleus collisions, the parallel ensemble method is commonly used, since it is 
simpler and numerically less expensive (see discussion in Sect. Ill C below). The advantage of the parallel ensemble 
method is that each parallel ensemble can be considered as a physical event, but this method does not converge to the 
solution of the kinetic equation in the limit of large N. 

This article is structured in the following way. First we introduce our GiBUU |33j ] transport model emphasizing 
the most relevant issues. Next, we discuss the solution of the Boltzmann equation using both the full ensemble and 
the parallel ensemble scheme and check their consistency. Finally, we present our results on DCX in comparison to 
the data. 



II. THE GIBUU TRANSPORT MODEL 



Boltzmann-Uchling-Uhlenbcck (BUU) transport models are based on the Boltzmann equation, which was modified 
by Nordheim, Uehling and Uhlenbeck to incorporate quantum statistics. A brief description of the formalism is given 
below. For a detailed discussion concerning the physical input for pion-induced reactions we refer the reader to Q 
and the references therein. A more general description of the whole GiBUU model will be given in a forthcoming 
paper |34j . 



A. The BUU equation 



The BUU equation actually consists of a series of coupled differential equations, which describe the time evolution of 
the single-particle phase-space densities f a (f,p,t). The index a = tt, cj, N, A, . . . denotes the different particle species 
in our model. A large number of mesonic and baryonic states is actually included, but at the energies of interest for 
this study, the relevant ones are w,N and the A(1232) resonance. 

For a particle of species X, its time evolution is given by 

dfx{r,p,t) _ dfx (r, p, t) | dH x dfx (f, p, t) dH x dfx (f, P, t) 
dt dt dp df df dp 

= Icoll (fx, fa, fb, ■ ■ ■) (1) 

with the one-body Hamilton function 



Hx(f,p) = \j(p + A x (r, p)) 2 +m\ + U x (f, P) + A° x (f, p) . 

The scalar potential Ux and the vector potential A x of species X may in principle depend upon the phase space 
densities of all other species. Hence, the differential equations are already coupled through the mean fields. In the 
limit of Icon = 0, cq. Q becomes the well-known Vlasov equation. The collision term I co u on the right-hand side 
of eq. Q incorporates explicitly all scattering processes among the particles. The reaction probabilities used in 
this collision term are chosen to match the elementary collisions among the particles in vacuum. Within the BUU 
framework the TtN reaction cross section is given by an incoherent sum of resonance contributions and a direct, i.e. 
point-like, contribution. Interference effects are therefore neglected. 



3 

B. Elementary processes 

DCX emerges from the interplay of two elementary mechanisms: pion-nucleon quasielastic scattering, with or 
without charge exchange, and pion absorption inside the nucleus. In our model, the cross section for quasi-elastic 
scattering is given by an incoherent sum of background a BG and resonance contributions 

0ViV->7riV = &7rN — >R — >7tN ®^rrN — >7rN ' 

The resonance cross sections are obtained from the partial wave analysis of ref. [36[. The background cross sections 
denoted by a BG are chosen in such a manner, that the elementary cross section data in the vacuum are reproduced. 
Background contributions are instantaneous in space-time, whereas the resonances propagate along their classical 
trajectories until they decay or interact with one or two nucleons in the medium. 

As an improvement in the earlier treatment 0, I^tI l3ftl | , we now included a more realistic angular distribution for the 
clastic scattering of the pions. Due to the P-wave nature of the A(1232) resonance, we assume for irN — > A — > nN 
in the resonance rest frame a distribution of the pion scattering angle 9 according to 

/ A M) = (1 + 3 cos 2 (8)) g(s, 9) 

which is peaked in forward and backward scattering angles. The function g(s,9), depending on Mandelstam s, 
parameterizes the energy dependence of the TtN angular distribution. In a coherent calculation the angular distribution 
is generated by interference effects, which can not be accomplished by our transport model. In our ansatz we need to 
split the cross section in an incoherent way to preserve our semi classical resonance picture. Therefore we take 

g(s,6) = {a- cos(9)f( mA -^ /m * 
with the A pole mass toa = 1.232 GeV. For the background events we assume 

r G ( Sl 8)=g( Sl 9). 

The constants a = 1.9 and /3 = 26.5 are fitted to the angular distributions measured by Crystal Ball [s^; a comparison 
of our parameterization to this data is shown in fig. ^ 
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FIG. 1: The angular distributions for the charge exchange process 7r~p — > ir°n in the CM-frame of pion and proton. The plots 
are labeled by the kinetic energies of the pions in the laboratory frame. The data are taken from ref. 




4 



P%fm~ 3 } 


o r t — 3i 
Pn[fm \ 


R p [fm] 


Rn[fm] 


a p [fm] a n [fm] 


0.0708 


0.0835 


5.194 


5.358 


0.4743 0.4780 



TABLE I: Parameters of the Woods-Saxon parametrizations for 45 Rh. 



Concerning pion absorption, the most important mechanisms are a two step processes in which irN — » A is followed 
by AA" — > NN or ANN — > NNN, and a one step background process irNN — > AW. As medium modifications, 
we include Pauli-blocking and Fermi-motion of the nucleons, Coulomb forces and hadronic potentials for baryons; we 
account as well for the collisional broadening in the A resonance width. In Q we discussed the influence of the real 
part of the pion self energy on absorption processes in nuclei. Since DCX is mostly sensitive to the surface (see the 
discussion in sec. II V|) . the influence of this real part has turned out to be negligible due to the low effective density 
and has not been included in the calculations presented here. For details see ref. 0] and the references therein. 



C. Nuclear densities 



The nucleons are initialized in a local density approximation. For g 6 0, ^Ca, 82 8 -^ we have implemented density 
profiles p(r) according to the parameterizations collected in ref. [39l |. which are of Woods-Saxon type for heavier nuclei 
{Ca^Pb) and of harmonic-oscillator type for lighter (O) ones. The proton densities are based on the compilation of 
ref. [401 from electron scattering. The neutron densities are provided by Hartree-Fock calculations. For \^ 3 Rh, we use 
a Woods-Saxon density distribution 

P° 

Pn(r) = 



1 + exp((r - Rn)/a„) 
P° P 

Pp(r) = l + exp((r P -R p )/a p ) ' (2) 

with the parameters given in table I |41|. The larger neutron radii of heavy nuclei play a relevant role in DCX as we 
show in Sec. IIVI 



III. SOLUTION OF THE BOLTZMANN EQUATION 



The fact that the DCX reaction depends considerably on the spatial distributions of protons and neutrons implies 
that it is also sensitive to the degree of locality of the scattering processes. In the non-discretized version of the 
BUU equation, the interactions are strictly local in space-time. Utilizing the so called test-particle ansatz to solve the 
problem numerically, this is no longer the case. Therefore, we elaborate in this section on this degree of locality of the 
scattering processes in our simulation. As a first step, we point out the connection of the underlying BUU equation 
to the actual numerical implementation. This will lead us to a proper definition of a typical volume in which particles 
are allowed to interact. Hereafter, we will introduce an approximation to the full solution: The parallel ensemble 
scheme. In this scheme computation time is radically decreased b y re ducing the locality. It is widely employed to 
solve the BUU equation in particle physics (e. g. |2(j, 0, [2^, I23I l37j) since otherwise computations would not be 
feasible. In sec. IIVI wc will discuss a possible problem related with this scheme and evaluate its applicability to the 
DCX process. 



A. The collision term 



For the sake of simplicity, let us consider a model with only one fcrmionic particle species which has no degeneracy, 
allowing only for binary scattering processes 

A(p~a) B(pb) — ► a(p a ) b(p b ) . 

If such a scattering event occurs at point f, then the single-particle phase-space density decreases in the vicinity of 
phase space points (f,p\) and (f,ps) and increases in the vicinity of (r,p^) and (r,p~b). Therefore, at each phase 
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space point, the collision term in eq. Q consists of a gain term due to particles which are scattered into this phase 
space point and a loss term due to particles which are scattered out: 

dmPXt) = I s ^(f,P%t)-I lom (r,PA,t) 



with 



dt 



I\ OSB (r,PA,t) = — ^ / d 3 p B [ dn C M d<Jj ^^ ab VAB f(f,p~A,t) f{r,ph,t)V a V b , 
{^) J J dilcM 



/gain {r, P% t) = ^—3 J d p a —-^ J d Pb{2ir) Vab f(r,Pa,t) f{r,Pb,t)VAPB 



where da AB^ab/^CM and do a b^AB / d?PA are the angular and momentum differential cross sections for the reactions 
AB — > ab and ab — > AB, respectively, vab and v a b are the relative velocities of the collision partners AB and ab. The 
terms Vx — 1 — f(r,p~x, t) with X = A, B,a,b correspond to the Pauli blocking of the final states. The momenta of 
the particles in 7i oss and /gain satisfy the condition of energy and momentum conservation. 

To point out the connection between our numerical implementation and the underlying BUU equation, we concen- 
trate on the loss term of BUU. We will therefore not elaborate on the gain term / ga in which describes the production 
of particles. However, its numerical implementation is analogous to the loss term since both are connected by detailed 
balance. 

The Boltzmann equation can be solved numerically using the so called test-particle ansatz where the single-particle 
density is expressed in terms of ^-functions 



f{r,p,t) = Jim 



N 



3 AxN 



]T 5(r-n(t))S(p-Pi(t))- 



Here A denotes the number of physical particles and N is the number of test-particles per physical one. So the single- 
particle phase-space density is interpreted as a sum of all test particle densities. The centroids of the ^-functions fi and 
Pi obey the classical Hamiltonian equations. The change of the single-particle phase-space density per infinitesimal 
At due to collisions is given by 



Af(r,p,t) = Ai(/ gain -/iosb) 
In terms of the test-particle ansatz the loss term reads 



(3) 



At Ji oss (f, Pa, t) 



lim 

N— >oo 



(2tt) 

N 



3 AxN AxN 



V V S(pA-pMr-fi) lim — dncMPaPb-^f 
r-i jv->oo aij J dilc 



ij — >ab 



M 



i = i 

3 ^ i 



=v a v b 



XvijAtVij—Sir 



= lim 



t 

(27T) 



N 



i) 



3 AxN AxN 



V V S(p A - Pi)S(f- fifPM / 5{r' - r 3 )d\' 

1=1 . , JAV 'i 



\ 



3 = 1 
3 + * 



I 



(4) 



(5) 



where and Vij are the total interaction cross section and the relative velocity of the test-particles AVij = 
UijAtVij/N is an infinitesimal volume in the vicinity of r^. Note, that the latter volume defines the locality of the 
scattering process of two test-particles. The term V a Vb denotes the blocking of the final state averaged over its angular 
distribution. We excluded self-interactions - therefore a test-particle cannot scatter with itself. 
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B. Numerical implementation 



In a real calculation the number of test-particles TV is chosen finite, for our purposes usually of the order of 300-1500. 
The time step At is chosen such that the average distance travelled by the particles during At is less than their mean 
free path. Therefore, AV^- is small enough such that a particle has no more than one scattering partner at a given 
time step. The algorithm proceeds as a sequence of the following steps: 

• First, we propagate at each time step the test-particles according to Hamilton's equations. 

• The loss term is implemented according to equation JSJ. Therefore we consider each term in eq. (JSJ separately. 
For simplicity let us just consider one summand describing the loss of the zth test-particle due to a collision 
with the jth. 

1. The term J AV .. 5(r ' — rj)d 3 r' gives 1 or depending on the fact whether j is within AVy. The volume 
AVij is chosen to be a cylinder of height Atvij with a circle basis Uij/TV; the symmetry axis is chosen along 
Vij and the basis is centered at fj . This corresponds to the usual minimum distance concept ■ 

2. If the result of the integral is 1, then we evaluate 



1 f ,„ da 



V a Vb = — / dilcM VaV b 



ij — >ab 



dn 



CM 



For this we make a Monte-Carlo integration with only one integration point, which is a good approximation 
in the large TV limit. This one point f2 CM is chosen in the center of mass (CM) frame randomly according 
to the weight --r~ S'"" 11 . Since y/s is fixed, f2 CM defines the random momentum p a CM . Furthermore, 



Finally, by boosting the momenta to the computational frame, we get 



V a V b = (1 - f(f,p a ,t))(l - f(r,p b ,t)) 



3. Now we interpret V a Vb as a probability that the reaction takes place. So we make a second Monte-Carlo 
decision on whether we accept the reaction or not. This corresponds to substituting V a Vb by a Bernoulli 
distributed random number with p = V a Vb- With this substitution, the expectation value of the summand 
is equal to the original summand. In the limit of many ensembles TV, i.e. many summands, this yields the 
right loss term. 

4. If the reaction is accepted, then wc get for this event where i is scattering with j the loss contribution 

5{p A -Pi)5(F- n) 

which corresponds to the destruction of the ith test-particle. Due to the double sum in eq. (JSJ, we get also 
the contribution 



S(pA-Pj)S(f- fj) . 

This latter term corresponds to the destruction of the jth test-particle. Note that we do not evaluate 
V a Vb for this case, but take the same value which lead to the destruction of i. This reflects that energy is 
conserved on an event-by-event basis. 

5. In our simulation, the final states with momenta pb CM and p a CM contribute to the gain term. New test 
particles with those momenta are therefore added to the simulation. 

A generalization to 2 — > 3 and 3 — > 2 processes and to finite particle species including degeneracies is straight- 
forward. 



C. Full and parallel ensemble method 



The kind of simulation wc described in the last section, is called a full ensemble calculation. There exists a common 
simplification to this method: the parallel ensemble method |22j| . In this scheme one sets TV = 1, performs TV runs at 
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P° P [fm- 3 ] 


pi [fm- 3 ] 


Rp [fm] 


Rn [fm] 


a p = a„ [fm] 


default parameterization ^39] 


0.0631 


0.0859 


6.624 


6.890 


0.549 


old parameterization [2] 


0.0590 


0.0900 


6.826 


6.826 


0.476 



TABLE II: Parameters of the Woods-Saxon parametrizations for §2 8 Pb. 



the same time and then averages the results over all runs. The densities used in each run are the averaged densities 
of all ./V parallel runs. Therefore the propagation part stays the same, whereas the collision term gets very much 
simplified. 

Note that the only justification for this simplification is a great gain in computation time. In a full ensemble 
method, the propagation part scales according to the number of test-particles per nucleon N, whereas the collision 
term scales with N 2 - therefore the computation time is 0(N 2 ). In a parallel ensemble method TV runs are performed, 
which results in O(N) computation time. So there is linear scaling in a parallel ensemble run, but a quadratic one 
in a full ensemble run. In pioneering works, it was shown by Wclke et al. |3l| and Lang et al. |30| . that the parallel 
ensemble scheme is a good approximation to the full ensemble scheme under the conditions of high-energy heavy ion 
collisions. 



IV. RESULTS OF THE SIMULATIONS 
A. Comparison of full and parallel ensemble runs 

In the previous section we introduced the concept of the parallel ensemble approximation. For DCX, surface effects 
are expected to be important, therefore the spatial resolution could be relevant in this context. Indeed, a major 
problem of the parallel ensemble scheme is that the volume AUy can become very large. In the energy regime 
under consideration, the incoming pions interact strongly with the nucleons so that the total cross section can reach 
more than 200 mb. For a parallel ensemble run, the typical volume has therefore the size of 5 fm 3 for a typical 
At = 0.25 fm/c. Since it is not obvious that the parallel ensemble scheme should be reliable in this regime, we hence 
decided to evaluate this approximation scheme by comparison to the full ensemble method. 

In fig. |21 the results for da/dfl at 180 MeV kinetic energy of the pion are presented for both methods. The results 
obtained are consistent with each other, therefore we used the parallel ensemble method to save CPU time and to 
reach higher statistics for all further calculations. Note that in the present problem, in order to obtain a result at a 
given energy for one specific nucleus one CPU-day is required in the parallel ensemble scheme. In the full scheme this 
takes of the order of 30 CPU-days for an acceptable statistics, as shown in fig. 

B. Influence of the density profile 

The DCX is, due to the low pion mean free path in nuclear matter Q, very sensitive to the surface properties of the 
nuclei. Therefore, we compared the results with our present density parameterization for 208 Pb [3j|, as described in 
section lTl CI to the results obtained with the one used in previous publications QIUUll- For 208 Pb, both distributions 
are parametrized according to eq. with the parameters given in table II. However, neutron skins are very interesting 
because in those skins only tt + mesons can undergo charge exchange reactions. For the positive pions this causes an 
enhancement of DCX processes at the surface, so the pions do not need to penetrate deeply for this reaction. Hence, 
the probability for their absorption is reduced. As can be observed in fig. 3, the enhancement in the cross section 
due to the neutron skin is roughly 35% at 180 MeV. We conclude that surface effects are prominent and can not be 
neglected. Note that for the tt~ similar arguments lead to a reduction of the cross section. 

C. Comparison to data 

Now we proceed to the comparison with the data measured at LAMPF by Wood et al. [l^. We discuss first the 
total cross section. Hereafter, we will explore angular distributions and, finally, the double differential cross sections 
as a function of both angles and energies of the outgoing pions are addressed. 

In fig. 0] one can see the excellent quantitative agreement to the total cross section data at 120, 150 and 180 MeV 
for Oxygen and Calcium. For the Lead nucleus we see some discrepancies. Notice that we reproduce the different A 
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FIG. 2: Comparison of the full and parallel ensemble methods. Angular distributions for the double charge exchange process 
± A -> -n^X at Skin = 180 MeV. The error bars denote the la statistical error of each point of the full ensemble run; the thin 
lines denote the la statistical error of the parallel result. 



dependencies of both (7r + ,7r~) and (7r _ ,7r + ) reactions. It is due to the fact that, when A increases, the number of 
neutrons increases with respect to the number of protons, and this favors the 7r + induced reaction. 

In [l^, Gram et al. discuss a scaling law of the total cross section. They argue as follows. Since the first collision 
takes place predominantly at the surface, the cross section should scale with A 2 / 3 . Furthermore they assume that 
DCX is mainly a two-step process, and that a pion which undergoes an elastic process at the first collision will not 
contribute. This is reasonable because the incoming pions loose energy in the elastic process, and their cross section 
for a second charge-exchange reaction is hereafter very much reduced. For a negative pion the first charge exchange 
reaction occurs with a probability of Z/N where Z(N) denotes the number of protons (neutrons). This is the case if 
the interaction is dominated by the A resonance, as it should be in this energy region. Finally, the second charge 
exchange process then takes place with the probability (Z — 1)/(A — 1) since, in the isospin limit, the ir° interacts 
equally well with protons and neutrons. Putting these considerations together and extending them to the 7r + case, 
the cross section for DCX is expected to scale according to 

A 2/3 Q Q-i ff{ ^ 
atot ~ A A^Q—> (6) 

where Q denotes the number of protons in the case of tt~ induced and the number of neutrons in ir + induced DCX. 

Gram et al. |19| find good agreement of this scaling law with experimental data. Also in the GiBUU simulation this 
scaling is fulfilled as can be seen in fig. El Nevertheless, one may wonder why this scaling law works in a process which 
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FIG. 3: Influence of the density distribution on the angular distributions for the double charge exchange 
at -Ekin = 180 MeV. The solid line shows the result obtained with our present density distribution 1391 
obtained with the previously used one [H I2I l37| . which contains no neutron skin. 
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FIG. 4: The inclusive double charge exchange total cross section as function of the nuclear target mass at i?kin = 120, 150 
and 180 MeV. The lines connecting our results are meant to guide the eye; the data are taken from Ref. (left panel: 
E kin = 120, 150 and 180 MeV, right panel: only 180 MeV). 



is so sensitive to the neutron skin on heavy nuclei, as has been shown in fig. 3. Since the first collision takes place 
on the surface, a neutron skin causes an enhancement in the A(ir + ,tt~)X reaction while A(n~ ,n + )X is suppressed. 
This effect leads to a deviation from the scaling. However there are also Coulomb forces which are not negligible. 
The Coulomb force enhances A(tt~ ,tt + )X by attracting the negative projectiles and repelling the positive products, 
which therefore have a smaller path in the nucleus and undergo less absorption. And, due to similar arguments, the 
reaction A(tt + ,tt~)X is suppressed. We find that this effect counteracts the one from the neutron skin restoring the 
scaling. In any case, the approximate scaling exhibited by the cross section shows that the reaction is very much 
surface driven and can be very well understood in terms of a two-step process. 

In fig. El we show da/dfl for DCX at E kin = 120, 150 and 180 MeV on 16 0, 40 Ca and 208 Pb as a function of the 
scattering angle 9 in the laboratory frame. Our results (bold lines) are shown together with their uncertainties of 
statistical nature (thin lines). The latter ones are well under control except at very small and very large angles, where 
statistics is very scarce. Again, there is a very good quantitative agreement for both O and Ca. In the Pb case, the 
(tv~ , 7r + ) reaction is well described, but the (ir + , 7r~) one is underestimated in spite of the enhancement caused by the 
neutron skin. 

Going into further details of the energy distribution of the produced pions, we show in figures [7|[S] and the results 
for da / (dfl dE^in) at different laboratory angles 9, as a function of the kinetic energy of the outgoing pion -Ekin- The 
overall agreement is good, better at forward and transverse angles than at backward angles. We observe a lack of 
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A -> X Jt" /+ 

: 1 1 1 r~- 

ji + induced + 
it" induced O 

; 9 ; 

! • O 

i i i i 

8 20 82 126 

A-Q 

FIG. 5: The scaling of the total charge exchange cross section according to eq. JJjJ is visualized by dividing atot by the factor 
A 2 ' 3 Q(Q — 1)/{A — 1) and plotting it as a function of A — Q. Q denotes the number of protons in the case of (tt~ , ty + ) and the 
number of neutrons in (7r + ,7r~). The points are GiBUU results at pion kinetic energies of 180 MeV; the dashed line denotes a 
function proportional 1/(A — Q), corresponding to the exact scaling. 

pions with energies below Ekin— 30 MeV. This feature becomes more prominent when going from O to Pb and is 
present for both incoming tt + and tt~ . The same problem also shows up in the work of Vicente et al. [TTj | (see their fig. 
9). A solution to this problem is not clear. Due to the low-energy nature of those missing pions, one may speculate 
whether quantum mechanical effects are responsible for the enhancement, and therefore can not be described by a 
semiclassical transport theory. 

V. SUMMARY 

We have studied pionic double charge exchange on different nuclear targets ( le O, 40 Ca and 20S Pb) in the A region 
(-Ekin = 120, 150, 180 MeV) with a semiclassical couple channel transport model (GiBUU). 

We have established the validity of the parallel ensemble scheme for this reaction, which is very sensitive to local 
density distributions by contrasting the results with those obtained in the more precise but time consuming full 
ensemble method. 

Furthermore, we compared the results of our model with the extensive set of data taken at LAMPF 0|, achieving 
a good agreement, not only for the total cross section, but also for angular distributions and double differential cross 
sections. Still, we miss some strength at backward angles and pion energies below -Ekin ~ 30 MeV. The scaling of 
the total cross sections pointed out in [l9l | could be reproduced. However, we found that two important effects that 
break this scaling: neutron skins and Coulomb forces compensate each other. 

We conclude that the implementation of pion rescattering and absorption in the GiBUU transport model successfully 
passes the demanding test of describing double charge exchange reactions. Thus the semi-classical approach is well 
suited to describe pion dynamics in nuclei for pion kinetic energies greater Ekin ~ 30 MeV. 
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FIG. 6: Angular distributions for the double charge exchange process n ± A — > n^X at E^i n = 120, 150 and 180 MeV. The data 
points are taken from [lr3 |: only systematical errors are shown. The solid lines represents the GiBUU results at 120 MeV, the 
dashed lines at 150 MeV and the dotted lines at 180 MeV. The bold lines represent the results while the thin lines represent a 
la confidence level for each point based upon our statistics. 
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